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ABSTRACT 


Presented  in  this  report  are  finite-difference  heat-transfer  equa¬ 
tions  for  transient,  radial  heat  flow  in  spheres  and  cylinders  and  for 
transient,  one-dimensional  heat  flow  in  flat  plates.  The  derived  equations 
apply  to  structures  before,  during,  and  after  surface  recession  for  all 
three  basic  structure  configurations  and  for  several  generic  material 
skin  combinations. 

For  each  skin  configuration,  the  accuracy  of  the  finite-difference 
procedure,  compared  with  exact  analytical  methods,  depends  on  opti¬ 
mum  selection  of  the  calculation  time  increment  and  the  incremental 
distance  between  temperature  nodes  in  relation  to  the  material  thermal'' 
properties  and  on  the  closeness  of  the  approximate  temperature  gradi¬ 
ents  to  the  true  gradients.  In  addition  to  these  common  criteria,  the 
magnitude  of  the  surface  recession  rate  in  relation  to  the  calculation 
time  increment  and  temperature  nodal  point  distance  affects  the 
accuracy  of  the  finite-difference  temperature  results.  When  compared 
with  exact  solutions  applicable  to  semi-infinite  flat  plates  undergoing 
surface  recession,  the  calculated  finite-difference  temperature  gradi¬ 
ents  during  recession  are  very  accurate  when  the  amount  of  material 
removed  during  n  calculation  time  increment  is  equal  to  or  less  than 
one  fourth  of  the  selected  distance  increment  between  temperature 
nodes. 

The  cylindrical  and  spherical  equations  are  presented  for  centrip¬ 
etal  heat  flow  and  surface  recession.  Two  simple  methods  of  converting 
the  centripetal  equations  to  the  centrifugal  form  for  applications  to 
structures  such  as  blast  tubes,  rocket  motor  combustion  chambers, 
and  nozzles  are  discussed.  These  two  methods  involve  making  a  minor 
number  of  sign  changes  in  the  centripetal  heat-flow  equations. 

Attractive  features  of  the  ablation- conduction  method  described  in 
this  report  are  the  negligible  increase  in  required  computer  time  over 
a  nonreceding  case  when  all  other  parameters  are  identical.  Secondly, 
the  nonshifting  temperature  grid  prevents  confusion  in  interpreting 
computer  results  and  readily  lends  itself  to  automatic  plotting  techniques. 
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SYMBOLS 


Area. 

Angle  in  radians. 

"olid  angle  in  steradians. 


Dimensionless  modulus 


R 

Z>n 


-  Radius  of  cylindrical  or  spherical  section. 


Distance  from  outer  surface  of  cylindrical  or  spherical 
section  to  temperature  point  n  {includes  thickness  of  ablation 
into  t  layer,  a). 


k 

c 

P 

T 

^net 

a 

L 

a 

Z 

B 
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Time  increment  for  computation. 

Temperature. 

Thermal  conductivity  of  material. 

Specific  heat  of  material. 

Density  of  material. 

Incremental  thickness  for  each  material. 

Net  heat  flux  at  boundary. 

Summation  of  ablation  into  any  one  t  {0  ^  a  <  t),  y^(aAt). 
Length  of  cylinder  (unity). 

Ablation  or  recession  rate. 

R  -]T)t2  . 

Symbol  notation  defined  after  each  use. 

Material  thickness. 


Subscripts  and  Superscripts 


a 

b 

c 

bs 

o 


-  Material  "A.  " 

-  Material  "B.  " 

-  Material  "C.  " 

-  Backside  or  internal  surface. 

-  External  surface. 
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-  Internal  surface. 

-  Conditions  existing  after  the  lapse  of  one  At. 

h  -  Nodal  point. 

m  -  Melt  condition. 


vm 


S*ctioni.  INTRODUCTION 


With  the  continuing  use  of  rel'able  ablating  reinforced  plastics  and 
subliming  materials  for  efficient,  economical,  thermal  protection  of 
missile  airframes  and  components,  an  accurate  simple  solution,  to 
problems  of  transient  heat  flow  in  solids  experiencing  a  variable  sur¬ 
face  recession  rate  at  one  surface  is  required.  This  solution  is  not 
only  necessary  for  flat  plates  but  also  for  cylinders  and  spheres. 

General  analytical  solutions  for  structures  undergoing  surface  reces¬ 
sion  are  not  available,  and  exact  solutions  are  known  only  for  special 
flat-plate  cases. 

The  analysis  of  small  hemispherically  tipped  vehicles  can  be  more 
accurately  calculated  by  a  spherical  program  than  a  flat  plate.  Often 
many  small  semicylindrical  leading  edges,  blast  tubes,  motor  cases, 
and  nozzles  can  better  be  assessed  by  a  cylindrical  procedure  than  by 
a  flat-plate  procedure. 

A  large  majority  of  the  materials  used  for  thermal  protection  of 
supersonic  missiles  possess  a  very  low  thermal  diffusivity.  As  a 
result  one-dimensional  heat  flow  in  flat  plates  and  radial  heat  flow  in 
cylinders  and  spheres  are  sufficiently  accurate  even  though  the  heat 
input  usually  varies  along  the  exposed  surface. 

A  numerical  finite-difference  method  for  heat  flow  before,  during, 
and  after  surface  recession  on  flat  plates,  cylinders,  and  spheres  is 
described  in  this  report.  The  equations  derived  for  cylinders  and 
spheres  are  for  centripetal  surface  recession;  however,  two  simple 
methods  of  using  the  same  equations  for  centrifugal  surface  recession 
are  discussed. 

A  brief  comparison  of  calculated  temperature  distributions  with 
exact  results  is  discussed  for  special,  ablating  flat-plate  cases.  In 
addition  to  the  criteria  affecting  the  accuracy  of  finite- difference 
results  for  a  piate  with  no  recession,  the  accuracy  of  the  numerical 
calculations  for  surface  recession  depends  quite  heavily  on  the  judicious 
selection  of  the  incremental  node  thickness  and  calculation  time  incre¬ 
ment  in  terms  of  the  actual  surface  recession  rate. 

The  advantages  of  the  ablation- conduction  method  presented  in  this 
report  are  the  simplicity  of  its  formulation,  the  versatility  of  the 
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boundary  conditions  (variable  recession  rate,  numerous  material  com¬ 
binations,  and  variable  thermal  properties),*  and  the  short  computer 
time  required.  For  the  same  structural  arrangement  and  identical 
selections  of  variables  such  as  node  thickness  and  calculation  time 
increment,  a  recession  computation  requires  a  negligible  increase  in 
computer  time  over  the  nonrecession  case. 

/ 
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*Tomperaturo  dependent  approximations  for  specific  heat  and 
thermal  conductivity. 
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Section  II.  HEAT  CONDUCTION  WITHOUT  SURFACE  RECESSION 


'1.  Transient,  One*  Dimensional  Heat  Transfer  for  Flat  Plates 

One-dimensional,  flat- plate  heat  transfer  in  a  homogeneous 
material  may  be  determined  by  solving  heat  balance  equations  at  the 
exposed  surface,  unexposed  surface,  interior  nodes,, .  and  interfaces. 

The  forward  finite-difference  method  was  used  to  solve  the  heat  balance 
equations.  It  was  assumed  that  the  incremental  thickness  (t)  can  be 
selected  small  enough  to  give  accurate  temperature  gradients  between 
adjacent  nodes  and  that  the  incremental  time  (At)  is  small  enough  to 
neglect  any  effect  on  regions  more  than  one  t  from  the  node  in  ques¬ 
tion.  The  stability  criteria  for  the  forward  finite  difference  equations 
can  be  found  in  Report  No.  RS-TR-65-  1. 1 

a.  Thick  Material 

(1)  Exterior  Surface.  From  Figure  1  the  heat  balance 
at  tl  exposed  surface  is 

qnetQ  "  qcond  =  qstor  -d 

l  —  z  l  (1) 

where 


qnPi-  =  net  heat  received  per  unit  area 
(T  T  ) 

qcond  =  *a  A 

'  o 


1  —  2 


„  „  la  .  (n -  T.) 

^stored  “  Pa  ca  i  A 


The  area,  A,  is  uniform  for  one-dimensional,  flat-plate  heat  transfer. 
Rewriting  Equation  (1)  we  have 

„  tl  lT  t  1  _  „  c  la  (T.  •  Ti) 

qnetQ  "  T  lTi-T2)  -Paca  2  {2) 


1U.  S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama, 
SOLUTION  OF  TRANSIENT  HEAT  TRANSFER  PROBLEMS  FOR  FLAT 
PLATES,  CYLINDERS,  AND  SPHERES  BY  FINITE- DIFFERENCE 
METHODS  by  W.  G.  Burleson  and  R.  Eppes,  Jr.,  15  March  1965, 
Report  No.  RS-TR-65- 1  (Unclassified  Report)  AD  461  662. 
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k-a  /_  „  \  Ta  (Tbs  "  ^hs) 

~  (Tbs-i  -  Tbs)  -  qneti  =  Pa  ca  ^ - ^ - 


Rearrange  and  solve  for  Tbs  with 
ka  At 

0  =  - § - 

Pa  Pa  ca  Ta 


Tbs  =  Tbs  +  ^Pa  (Tbs-i  “  Tbs)  "  ‘Ineti 


Pa  ca  Ta 


b.  Thick- Thick  Material 

At  the  interface  between  material  "A"  and  "B"  (Tn> 
Figure  2),  the  energy  balance  is 

^cond  +  ^cond  ~  ^stored 
n-i—n  n+l— n  n 


ka  ,  >  ku  . 

—  (Tn_!  -  Tn)+—  (Tn+1  -  Tn) 


=  (pa  ca  ~2 


ra  Tb\  (Tn  "  Tn) 

F+PbcbTj  - 


Material 

A 


Material 

B 


_ ^  *0  ^  ^  ^  ^  ^  Tb  ^ _ 


Figure  2. 


Rearrange  and  solve  for  Tn 


2  ka  At 


Tn  -  Tn  +  (Tn-t  "  Tn)  ra  (pa  ca  ra  +  pb  cb  Tb) 

+  (T  .  T  )  _ 2kb  At _ 

[  n+i  n)  T,/p  c  T  +Ph  ch  tu| 


rb(Pa  ca  Ta  +  Pb  cb  Tb) 
c.  Thin- Thick  Material 


(12) 


At  the  exposed  surface  (Figure  3)  the  energy  balance  for 
the  thermally  thin- thermally  thick  interface  is 


^netQ  +  ^cond  “  ^stored 
2—1  l 


(13) 


or 


kb  /_  _  \  /  Tb\  (t!  -Ti) 

qneto+^  (T2  -  T!)=  ^Pa  ca  ra  +  Pbcb— J  — 


(H) 


Rearrange  and  solve  for  Tj 

^netQ  At 


Ti  =  Ti + 


Pa  ca  Ta  +  Pb  cb 


_ ^b  At  ( Tz  -  Ti) 

Tfc)  ^Pa  ca  Ta  +  Pb  <=b-yj  (15) 


Material 

Material 

A 

B 

Thin  Skin 

qnet„ 


1 

- 1 - 

,TI  |  T2, 

1 

1 

1 

,  |  T3  ( 

1 

1 

_ ( 

Ta 

Tb 

— 

2 

Figure  3. 


d.  Thick- Thin  Material 

The  energy  balance  at  the  backside  surface  (Figure  4), 
T^g,  may  be  written  as 

^cond  ”  ^neti  ~  ^stored 

bs-l-*bs  bs  (16 


^a  /  \  /  Ta  \  (^bs  ”  "^bs) 

~  (Tbs-i  -  Tbs)  -  qnet^  =  ^Pa  ca  “J"  +  Pb  cb  Tbj  ^ - 

Rearrange  and  solve  for  T^g 

■  ka  At  (T’qs-1  ~  Tbs) 


Tbs  =  ^bs  +  /  ra  \ 

Ta  (  Pa  ca  +  Pb  cb  Tb) 

^neti  At 

-  7 - - -r  •  <18> 

faC*T*'bVk) 

e.  Thick- Thin- Thick  Material 

At  the  interface  between  material  "A"  and  "C"  (Tn, 
Figure  5),  the  energy  balance  for  the  thermally  thin  material  "B"  is 

^cond  +  ^cond  =  ^stored 

n- 1— n  n+t— n  n  (191 


+  Pb  cb  Tb 


tjt  (Tn.,  -Tn)t^  (Tn+i  -  T„)  =  |pa  ca-S 

x  x  tc\  (T;  -  Tn)- 

+  Pb  cb  Tb  +  Pc  cc  2  J  At 
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Material 


■SSrtt- 


i  fcvi,  ?*?' ' 


■s-^V.'l^-yy  , '  ~  -JV 


Ife 


II 


/ 


Rearrange  and  solve  for  Tn 

ka  At  ^Tn_t  -  1  n) 


Tn  + 


Ta 


Ta^Pa  ca*2  +  Pb  cb  Tb  +  Pc  c' 

kC  At  (Tn+1  -  Tn) _ 

V 


('a'af 


(21) 


+  Pb  cb  Tb  +  Pc  cc  Y 


1 

r) 

2, 

Transient,  Radial  Heat  Transfer  for  Cylinders 

a.  Thick  Material 

S 

. 

(1)  Exterior  Surface.  Consider  a 
heated  as  shown  in  Figure  6.  From  the  energy 
surface 

<Un 

A-i  "  ^out  =  ^stored 

1—2  1 

where 

<Hn 

=  ^net0  ’  “  R9L 

(22) 


^out 
1  —  2 


(Ti  -T2)  ,  Az  =  ^R  -  0L 

(*■■?) 


U  _  T  ) 

^stored  =  Pa  ca  “  »  ^3  =  ^  ~  ^ 

l 


or 


R0L  qnetQ  -  (r 

=  (r 


Ja 

2 


;)  6L  Ta  (T>  ‘  T*) 

t)  eL(Pa  c*t) 


Ta)  (t;  -  tJ 
2  /  At 


(23) 


Let 


ka  At 

Pa  =  — - r 

a  Pa  ca  Ta2 


i " 


Let 


p  =  . 
Pa  ca  V 


rearrange;  and  solve  for  Tn 

/  _  T 


-£n 


Tn  =Tn  +  2  )  (T»-<  -  T") 

fR  _y\  .  If 

“  n  2  1/  \ 
+  M  R  .  2t  J  (Tn+l  '  Tn) 


(27) 


(3)  Backside  Surface.  The  energy  balance  at  the  back¬ 
side  surface  (Figure  6),  Tbs>  ?"nay  be  written  as 


9cond  “  ^neti  "  ^stored  '^'3 
>  bs 


bs-l-*bs 

where 


^cond  '^1  =  r  (Tbs-i  "  Tbs 

bs-i— -bs 

'net: 


(28) 


)(■&.*?) 


0L 


a2 


'neti  (R  'X>bs) 


0L 


Ta  (^bs  “  Tbs) 
^stored  “  Pa  ca  2 
bs 


At 


(■2>.*t) 


0L 


or 


0L 


T  ^Tbs  +  2  ^  ("Tbs-i  "  Tbs)  "  9L  ^netj  "23  Tb^ 

Ta  Ta'\  (Tbs  ~  Tbs) 

T[*sy»*T)— — 


=  0L  Pa  c 


(29) 


Let 


P,  ■ 


ka  At 

Pa  ca  V 


11 


I 


I 


rearrange;  and  solve  for  T^s 

R  -£Tbs  +  . 

^bs  “  ^bs  +  ^PaJ  ~  Ta  I  ^bs-i  "  ^bsj 


R  -£' 


bs  T  4 


^net. 


2  At 


R  -ETbs 


Pa  ca  Ta  \ R  -£rbs  +  If 


(30) 


b.  Thick- Thick  Material 


At  the  interface  between  material  "A"  and  "B"  (Tn, 
Figure  7),  the  energy  balance  is 


0cond  +  ^cond  “  ^stored  ^3 
n-l— n  .  n+i— n  n 


whe  re 


0L 


•Icond  Al  *  ^  (Tn-i  -  Tn)  (R  -ST«  +  t) 


(t„+i  -  Tn)  (r  -£t„  -  9 L 

(R  -Et"  +  t)  Pa  c*  T 


n-l—  n 

^cond  ^2 
n+i—  n 


hb 

Tb 


^stored  -^3  = 
ri 


(R  'ETn  -  t) 


«  -  Tn) 


At 


0L 


(31) 


(32) 


12 


Figure  7. 


Solving  for 


Tn  *  Tn  + 


2  ka  At  (r  -^Tn  +t)  (Tn-1  '  Tn) 

Ta(R  'Z/n  +  t)  Pa  Ca  Ta  +  (R  £/n  "  T  ) 

2  kb  At  ~y  Tn  -  |Tn+i  *  Tn) 
”^)Tn  +  ~4~J  Pa  ca  Ta  +  ^R  ~])j^Tn  "  ~4fjpb  cb  Tb 


Pb  cb  Tb 


1  (33) 


c.  Thin-Thick  Material 

At  the  exposed  surface  (Figure  8)  the  energy  balance  for 
the  thermally  thin- thermally  thick  interface  is 

%etQ  +  ^cond  Az  “  ^stored  "^3 


where 


^net0  Al  =  qnet0  R  0L 


^cond 


a2  ^  (tj-t.)  (r  -  Ta-^) 


GL 


mmmmmmmmmmmrnm 


d.  Thick- Thin  Material 

The  energy  balance  at  the  backside  surface  (Figure  9), 
Tbs'  may  be  written  as 

qcond  Al  ’  ^neq  Az  "  ^stored  As 
bs-i-*bs  bs 


where 


9cond  Ai  =  ”  (Tbs-i  ”  Tbs)  (R  ‘T^bs  +  2  )°L 
bs-i-*bs  11 

'■neti  Aj  =  %ct|  (R  -ETbS  -  Tb)  eL 

^stored  A’  =  Pa  'af  (R  £/bs  +  t) 
bs  L 

/  XT'  Tb\  (Tbs  ~  Tbs) 
+  Pb  cb  Tb  ^R  V  /bs  i  I  At 


^stored  Aj  pa  ca  2  \R 
bs  L 


,.Mot«rlol 

A 


a  t: 


»  Moterlol 
B 

Thin  Skin 


9 


_1 - 


Figure  9. 


or 


e.  Thick-Thin-Thick  Material 

At  the  interface  between  material  "A"  and  "C"  (Tn, 
Figure  10),  the  energy  balance  for  the  thermally  thin  material  "B"  is 

^cond  +  ^cond  =  ^stored 

n- 1— n  n+i— *-n  n  (40) 


where 

^cond  -Ai 
n-  l— n 

^cond  ^2 
n+i-~ n 
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3.  Transient,  Radial  Haat  Transfer  (or  Spharas 


a.  Thick  Material 


’•  ’  (1)  External  Surface.  Consider  a  spherical  segment  heated 

as  shown  in  Figure  1 1.  From  the  energy  balance  at  the  peripheral  surface 

^net  A*  "  ^cond  Az  =  ^stored  Aj 

°  1-2  1  (43) 


where 

clnet0  Al  =  qnetD  R  ^ 


^cond  Az  = 
1  —  2 

^stored  Aa 
l 


4> 


(44) 


(The  average  area  terms  are  derived  in  detail  in  Report  No.  RS-  TR-65-  1 .  )2 


2Burleson  and  Eppes,  loc.  cit. 
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Let 


Pa 


ka  Afc 

Pa  ca  Ta2 


f 


rearrange;  and  solve  for  Ti 


(2)  Interior  Node.  The  energy  balance  at  any  interior 
point  n  of  a  homogeneous  wall  (Figure  11)  may  be  written 

Icond  -Ai  +  qcond  -^2  =  ^stored  -^3 

n-i~*-n  n+i— -n  n  (46) 


qcond  Ax  =  ~  (Tn.,  -  Tn)  (r  £rn  +  ~ 

n-i— n  a  L' 

n+i— -n  L 


(Tn  -  Tn)  [7  ^  \2  Ta2l 
^stored  A-3  =  Pa  ca  Ta  ^  ^_^Tnj  +  12  ^ 

n  L  J 


*[(*  ■Lt"  +  t)!  +  1t]  ~ 

+  *  [(r  -  £)*  ♦  T4]  £  (t.«  -  T„) 
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Let 


Pa  = 


ka  At 


a  "  Pa  ca  Ta2 
rearrange;  and  solve  for  Tn 


T  =  T 
1  n  x  n 


+  Pa 


(48) 


(3)  Backside  Surface.  The  energy  balance  at  the  backside 
surface  (Figure  11),  may  be  written  as 


^cond  Al  ”  ^neti  Az  “  ^stored  A* 

bs 


bs-i-~ bs 
where 


(49) 


9cond  Al  '  (Tbs-i  ”  Tbs)  ^  ^Tbs  +  2  )  +  1 


bs-i—  bs 


%et:  A* 


a 

12 


4> 


^netj 


(R  -ETb.) 


4> 


^stored  Aj  "  Pa  ca  £ 
bs 


Ta  ('•Pbs  ~  ^bs) 


At 


a 

48' 


4> 


or 


(R  5>.  *  t)  1  T>  r  (TbS-l  -  Tbs)  -  1m.t,  ♦  (R  -Zru) 


4> 


f-E— t) 


Pa  ca  — 


(^bs  ~  ^bs) 


At 


(50) 


Let 


ka  At 

Pa  OCT2 

Pa  ca  Ta 
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1 


I- 


wm*<Mt**  ’■•asMtertttfiiSIlte’ wtos*':*  - 


rearrange;  and  solve  for  Tbs 


^bs-i  "  Tbs) 


(51) 


At  the  interface  between  material  "A"  and  "B"  (Tn, 
Figure  12),  the  energy  balance  is 


^cond  +  ^cond  =  ^stored  ^3 

n-l— -n  n+l— -n  n 


(52) 


where 


z  xa2 


*condA‘  =^(Tn-l-Tn )  (R  +  t)  '  12 

kb  /_  _  \[/L  Ib\2  A^b! 

"2-r n " 2 )  12 


^cond 


n+l— n 
^stored  =  ^ i 


xb  (Tn+i  "  Tnj 


4> 


n 


(’ 

'(•E’.*t M 


a 

48 


Pa  ca 


or 


k_  , 

^i(T 

n-,  *  Tn)  f  « 

f'- 

/  48 

Ta 

Pa  ca  —  ♦ 

1  n+i 


(t„  -  Tn) 


Tb 

2  (  At 


(53) 
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Z1"-'  ^T> 

I*  ;Tn  J 


it  i 


/!' 
\\  , 


Figure  12. 


7" 

e 


i 

Rearrange  and  solve  for  Tn 

T.  T  _ 

(T- 

’  ^  lit*  7 


c.  Thin-Thick  Material 

At  the  exposed  surface  (Figure  13)  the  energy  balance 
for  the  thermally  thin- thermally  thick  interface  (Ti)  is 


C*net^  “^1  +  ^cond  ~  ^stored 


Pn  '  a  Tn 


TbV  Tb* 


Thick-Thin  Material 

The  energy  balance  at  the  backside  surface  (Figure  14), 
T^s,  riiay  be  written  as 

^cond  Al  "  ^net-  A2  "  ^stored  A3 

bs-i— bs  bs  /c;r\ 


where 


bsf.-bt  %a  (Tbs-1  '  Tbs)  (R  STbs  +  T )  +  -fz  * 

^neti  A2  =  ^netj  |R  "]£Tbs  “  Tbj2  4> 

^stored  A3  =|pa  c*~  |R  ~^)Tbs  +  "4" j  +  4I 

|/„  Tb\ 2  Tb2]|  (Tbs  -  Tbs) 

+  Pb%Tb  (R-^bs-Tj  +12  f - At 


Tb2  {  \Tbs  _  Tb< 


Material 

A 


!\/F 

1  \  I  I  i 


J  /C  /  ' 

I  t  ill 

Jim 

i*-To>  Tb 

-*—**»—* 


Material^" 

B 

^ThlnSkln 


Figure  14. 


or 


4  ^  |(R  'X>  4  r)  4  IT  |  (Tbs- ,  -  -  *  ‘-net,  ( K  -2  *0. 


♦  Pa  ca“ 


IS  4  t)' 


+  Pb  cb  Tb  R 


ETb\  Tb2 
Tbs  "  2  /  +  12 


21)  (Tbs  -  Tbs) 
At 


Rearrange  and  solve  for  TjjS 

ka  ^  (R  -J}bs  ♦  ^  +  Yi\  (Tb-«  ‘  Tbs) 

Tbs  lb* f  -|P;iC^[(R  j>bs , Pb:;;;|(R *  %\\ 

/  V  < 

Wt,  At  (r  *E/bs  -  TbsjJ 

p« r  (R  •  t)1  *  is  •  °b Cb Tb  j(R -Et»s t)!;S 


*  .  Tb* 
'  11 


e.  Thick- Thin- Thick  Material 

At  the  interface  between  "A"  and  "C"..(Tn,  Figure  15), 
the  energy  balance  for  the  thermally  thin  material  "B"  is 

^cond  '^1  *  Rcond  ’  ^stored 


n-i— ■  n  n+j  —  n 


where 


<Uond  ■'>  (V.  •  ■„)  («  JV.  y)  •  f: 


^  oik!  t’  (^n-  !  ‘  ^  n)  |(k  y^]Tn  ’  Tb  '  ~£ )  *  'jj  ^ 

■a « l  »n  *  * 

"M-r.d  jf.  y  \{H  E-  -  t.t  • 

*  — |(«E-  tf’SI  o.;|(«En 

M  *‘,ll  f1'  1 J 

I  *  4 •*  ft  ~\l 


Figure  15. 


~  9 


-  pb  %  Tb  !(r  -T/h 
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+  P, 


Section  III.  HEAT  CONDUCTION  DURING  ABLATION 


To  calculate  the  heat  transfer  in  a  material  when  the  exposed 
surface  of  the  material  is  heated,  an  energy  balance  must  be  performed. 
For  this  analysis  the  energy  is  considered  to  be  either  radiated  away 
from  the  exposed  surface,  conducted  into  the  cooler  interior  of  the 
structure,  or  stored  in  the  material  near  the  exposed  surface.  If 
sufficient  energy  is  stored  at  the  surface,  the  surface  temperature  will 
eventually  reach  a  critical  value.  This  value  is  usually  known  as  an 
ablating,  melting,  or  subliming  temperature  (Tm).  In  this  report  it 
is  assumed  that  the  ablating  temperature  Tm  is  known  or  calculable, 
and  this  temperature  remains  constant  while  ablation  is  in  process. 
Another  basic  parameter  required  once  the  exposed  surface  has  reached 
the  ablation  temperature,  Tm,  is  the  recession  rate  (ablation  rate)  or 
the  rate  of  material  removal.  It  is  assumed  that  the  ablation  rate  (a) 
is  known  or  can  be  calculated  for  any  given  increment  of  time,  but  may 
change  as  a  function  of  time. 

It  is  also  assumed  that  once  the  exposed  surface  reaches  the  melt 
temperature  (Tm)»  the  recession  rate  governs  the  amount  of  material 
removed.  The  material  properties  of  specific  heat  (C)  and  thermal 
conductivity  (k)  may  all  be  a  function  of  temperature.*  The  pyrolysis 
of  the  material  leaving  the  heated  surface  of  the  slab  has  been  left  out 
intentionally  because  of  the  complexity  of  the  problem.  However,  this 
parameter  can  be  included  in  the  heat  balance,  if  desired. 

Figure  16  shows  the  temperature  grid  arrangement  for  a  slab 
undergoing  surface  recession.  Nodal  points  Tz  through  Tn  do  not 
change  or  shift  positions  while  ablation  is  occurring  between  the  orig¬ 
inal  nodal  points  Ti  and  T2.  When  the  ablation  front  or  receding  front 
reaches  or  passes  Tz,  the  T2  nodal  point  temperature  assumes  the 
value  of  Tm>  and  T3  is  calculated  in  the  same  manner  as  was  T2  when 
the  receding  surface  was  between  nodal  points  Ti  and  Tz.  Once  the 
ablation  front  reaches  a  nodal  point,  the  ablation  distance  "a"  is  reduced 
by  the  value  of  t3  (0  <  a  <  Ta),  and  the  process  starts  over  again. 

Let  us  first  assume  that  Ti  and  T{  take  on  a  constant  value  Tm 
for  all  times  during  ablation.  The  primed  or  future  values  for  all 
T's  (T3  through  Ta  in  Figure  16)  can  be  obtained  by  ordinary  forward 
finite-difference  methods  since  each  of  these  nodal  points  is  a  full  Ta 


^Temperature  dependent  approximations  for  specific  heat  and 
thermal  conductivity. 


from  the  two  adjacent  nodal  points.  Thus  only  the  future  temperature 
at  one  nodal  point  (Tj)  is  not  defined.  Examination  of  the  geometry 
for  the  ablation  case  shown  in  Figure  16  in  light  of  the  pure  conduction 
condition  described  in  Figure  1  reveals  that  two  basic  conditions  must 
be  considered  for  the  ablation- conduction  case.  One  case  is  for 


_  ^  Ta  ,  .  ,  .  ,  Ta 

0  Sai  —  i  and  the  other  is  for  —  <  a  <  ra. 


The  energy  balance  for  the  T2  nodal  point  using  forward  finite- 
difference  approximations  for  a  slab  is  as  follows  when  "a"  is  equal 
T- 

to  or  less  than  —  . 


80 


Note: 


T| 1  T|  =Tm  during  oblotion. 
Figure  16. 


Recession  condition  Os  a<  — 


^in^  '  ^out  ^  "^stored  ^ 
1—2  2—3  2 

where  A  is  unit  area,  or 


Ta  -  a 


(Tm  -  T‘)  +  ^  (Tl  -  T*)  ^  •>»  =a  ra 


(64) 


(65) 


29 


Let 


ka  At 

Pa  =  7  ~  r~2 
Pa  ca  Ta 

rearrange;  and  solve 


for  T2 


Tz  =  t2  +  Pa  (T~'a")  (Tm  "  Tz)  +  Pa  (T3  '  Tz)  •  (66) 

Ta 

This  equation  is  valid  for  0  £  a  <  —  .  When  "a"  =  0,  Equation  (66) 

reduces  to  an  interior  node  equation  having  the  form  of  Equation  (6) 

Ta 

for  the  pure  conduction  condition  with  Tm  =  Tj.  When  "a"  =  —  , 
Equation  (66)  is  unstable  for  (3  >  V3  . 

Ta 

For  the  case  where  —  <  a  <  Ta ,  the  reduced  storage  associated 

with  nodal  point  2  must  be  considered.  The  energy  balance  for  the 
T2  nodal  point  using  forward  finite  differences  becomes 


<*in  A  -  %ut  A  =  ^stored  A 
1—2  2  —  3  2 


(67) 


where  A  is  unit  area,  or 

ka  1 


|Ta  -  aj  (im  "  Tss/  +  Ta 


Pa  ca 


(t3  -  t2) 

/3ra  \  (t^-  Tz) 
\  2  /  At 


Let 


ka  At 
Pa  ca  Ta2 


rearrange;  and  solve  for  T2 

„  „  4  (Tm  '  T,) 

T2  =  Tj  +  Pa  ~  rpr. 


' a)  n-  • a) 


(t3  -  t2) 


(68) 


(69) 
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'  i.'.  r<~%. 


Equation  (69)  is  discontinuous  or  unstable  at  "a"  =  ra,  regardless 
of  the  value  of  p.  This  instability  would  bring  about  problems  each 
time  "a"-*Ta  and  nodal  points  are  removed. 


To  eliminate  the  stability  problems  associated  with  Equations  (66) 
and  (69),  an  investigation  was  made  of  a  backward  finite-difference 
approximation  for  T2  when  ablation  is  in  process.  The  backward 
finite-difference  approximation  is  accomplished  by  priming  all  tem¬ 
perature  values  to  the  left  of  the  equal  sign  in  Equations  (65)  and  (68). 
This  is  possible  since  the  surface  temperature,  Tm,  is  known,  and 
all  other  temperatures  beyond  T2  are  easily  obtained  by  standard 
forward-difference  techniques.  Two  resulting  equations  in  terms  of 
T2  are  obtained  in  relation  to  the  position  of  the  receding  surface. 


or  letting 

k-,  At 

pa  = - - - 5 

PacaTa2 
and  solving  for  T2 


(72) 


or  letting 


ka  At 
Pa  caTaZ 


i 


<  a  < 


T 


a  • 


(73) 


It  can  readily  be  seen  that  Equations  (71)  and  (7  3)  are  continuous 
for  all  values  of  ''a"  (0  <  a  <  Ta)  and  are  a  marked  improvement  over 
Equations  (66)  and  (69). 

A  mid  finite-difference  approximation  was  also  investigated  for 
finding  T2.  This  mid-difference  method  is  a  better  approximation  to 
the  exact  solution  than  the  backward  finite-difference  approach.  How¬ 
ever,  the  additional  complexity  of  the  equations  and  computer  storage 
requirements  were  considered  unwarranted  for  the  additional  accuracy 
gained. 

Another  ablation- conduction  approach  that  has  been  used  success¬ 
fully  is  the  "shift"  method  shown  in  Figure  17.  An  interpolation 
routine  is  used  with  known  temperatures  Tj,  T2,  T3,  and  T4  to  get 
T2  and  T3  located  at  even  increments  of  Ta  from  T|.  The  general 
forward  finite-difference  interior  equation  is  used  to  get  Ti,  i.  e.  , 

Tj,  T£,  and  T3  are  used  to  get  T£.  Nodal  points  T3  through  Ta  are 
calculated  using  the  ordinary  forward  finite-differen*.e  equations  and 
original  temperature  node  locations.  With  the  prime  temperatures 
known  (Tj,  T£,  T3),  it  is  possible  to  use  a  three-point  interpolation 
routine  to  find  T2  located  at  a  distance  of  one  Ta  in  front  of  nodal  point 
T3  and  Ta  -  a  distance  from  the  receding  surface.  This  "shift" 
method  can  be  used  accurately  until  the  receding  surface  becomes 
closer  than  two  Ta’s  from  the  unheated  surface.  At  this  time  other 
special  equations  must  be  utilized. 


After  consideration  of  the  advantages  and  disadvantages  of  the 
three  approaches  discussed,  it  was  decided  to  pursue  the  backward 
finite-difference  method  in  deriving  the  energy  balance  equations  (T2) 
during  ablation  for  one-dimensional  flat  plates  and  radial  conduction 
in  cylinders  and  spheres. 

The  finite-difference  equations  presented  in  this  report  for  finding 
the  temperature  profile  near  the  receding  surface  are  limited  to  the 
condition  of  0  <  a  <  Ta.  However,  as  the  receding  front  passes  succes¬ 
sive  originally  selected  temperature  nodes,  the  equations  as  derived 
are  applicable  if  appropriate  temperature  subscripts  are  used.  For 
example,  when  the  receding  front  is  between  the  original  location  of 
T2  and  T3  (Figure  16)  on  a  semi-infinite  slab,  either  Equation  (76)  or 
(79)  is  used  to  find  T3  by  increasing  all  temperature  subscripts  by 
one. 


1.  Flat  Plat* 


The  procedure  for  calculating  heat  flow  during  surface  reces 
sion  is  described  for  all  expected  material  combinations  with  the 
general  thick  equations  being  a  repeat  of  Equations  (70)  through  (74). 


a.  General  Thick  with  6a  >  Ta  (Figure  16) 


(1) 


0  <  a  < 


The  energy  balance  for  T2  is 


9in  "  9out  "  ^stored 
1  —  2  2  — *  3  2 


Letting 


Pa  = 


ka  At 
Pa  ca  Ta2 


(74) 


(75) 


and  solving  for  T2 


14 


; aa* 


Ta 

(2)  —  <  a  <  ra.  The  energy  balance  for  T2  is 


^in  "  ^out  ^stored 

1—2  2—3  2  (  (77) 


It  should  be  noted  in  Equation  (79) that,  as  "a"  approaches  Ta,  T2 
takes  on  the  value  of  Tm.  This  is  true  under  actual  conditions  since 
the  nodal  point  Tj  is  moving,  and  all  the  other  nodal  points  are  fixed. 

As  the  ablation  front  reaches  the  original  interior  nodal  points,  these 
nodes  take  on  a  temperature  value  of  Tm. 

As  long  as  the  remaining  wall  thickness  of  the  material  undergoing 
recession  is  equal  to  one  or  more  than  one  t,  the  interior  and  interface 
equations  presented  in  Section  II  are  used  to  solve  the  heat  balances 
throughout  the  structure  away  from  the  receding  surfaces. 


3o 


b.  Special  Thick  with  6a  s  Ta  (Receding  Surface  £  Distance 
Ta  from  Backside,  Figure  18) 


When  the  receding  surface  is  less  than  one  incremental 
ra  from  another  material  surface  or  the  backside  of  material  "A," 
special  consideration  must  be  made  for  all  boundaries  normally 
experienced. 


^  “a  ■ 

Figure  18. 


(1)  0  <  a  <  —  .  The  energy  balance  for  T2  is 


qcond  "  qneti  ~  ^stored 
l  —  z  2 


(80) 


or 


Letting 


It  t'\  „  o  ,  Ta  -  Ill 

\Tm  ‘  Tz/  "  qneti  -  Pa  ca~  ~ 


ka  At 

la  -  - 

a  Pa  ca  V 


(81) 


and  solving  for  T2, 


I 

T2 


/Ta  • a  \ 

^  ^neti  At  /Ta  -  a  \ 

\  Ta  / 

fz  ’  Pa  Ca  Ta  \  Ta  / 

2Pa  + 


(82) 


(2) 


The  energy  balance  for  T2  is 


-  <  a  <  ra. 


qcond  '  qnetj  ^stored 
1  —  2  2 


(83) 


or 


?7hr  (T-  -  T0  -  w,  -  pa 


(Ta  *  a) 

Letting 


2  -  T2 
At 


k_  At 

_  d  _ 

Pa  ~  0  c  T  2 
Pa  ca  Ta 


and  solving  for  T2  , 

^Pa  + 


Ta 


T2  - 


^  9neti 
Pa  ca  Ta 


At 


/Ta  -  a\  2 
2Pa  +  — 


(84) 


(85) 


c.  Special  Thick- Thin  with  6a  <  Ta  (Figure  19) 


m —  8a  — H 


Figure  19. 
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(1)  0  <  a  <  ~  .  The  energy  balance  for  T2  is 


or 


Pcond  "  ^neti  '  ^stored 

1—2  2 


aj  (Tm  ’  Tz)  '  qneti  :  (Pa  Ca  2 


(86) 


(87) 


Solving  for  T,  , 


ka  At 

•  *  T»  ('*  '  1 

At 

/Ta  -  -\ 

Taj 

Ta 

Pa  ■  a  y  '  Pb  <b  Tb 

1  1 

Ta 

Pa  Ca  y  •  P|)  %  t^ 

|  \  T*  ) 

*»  At 


TJ  -  ,1 


‘Pb‘bTb| 

Ta 

(2)  —  <  a  <  Ta.  The  energy  balance  is 


(88) 


9cond  "  ^netj  =  ^stored 

.  —  2  2 


(89) 


or 


/  »  \  /  T 

\  ni  "  V  "  qncti  "  (Pa  Ca  “ 
(tz  -  Tz) 


a  ■  a 


+  Pb  cb  Tb 


\  (t ii_ 
7 


(90) 


Solving  for  T2 


*  m 

„  !u  •  *\ 

'Imi,  At 

TJ 

Ta  -  a 

Pa  ca  ““2 —  *  Pb  *b 

rT,r-;  ■  i 

Ta  "  a 

Pa  ca  ^  ~  ♦  Pb  1  b  Tb 

At 

Ta  -  a 

Ta  | 

Ta  *  a  y 

Pa  ca  _»  *  Pb  cb  Tb  a 

(V) 


(91) 


(1)  0  «aS  —  .  The  energy  balance  is 


1 


'■■■I 


«  V  *■> 


^cond  +  ^cond  ^stored 


1  —  2  3—2 


TT^)  (Tm  ‘  t;)+  k  (Ti  ‘  Ti) s  h Ca^  ' Pb  Cb  Tb 


+  Pc  Cc 


Tcl  (Ti  -  T2) 


Solving  for  T i, 


/Ta  •  a\ 

\-^r) 


r  r<i 

l*'.  I  '  ..  £ 


pti  b  ti>  •  p. 


_ _ In^a*' _ _  T,  ke  At  _ 

[p.,  ■  ■,  y  Pb  b  Tb  -  p,  < ,  T  ]  4«'aT‘  P».  ‘  b  T>,  *  Pt  «c  t] 

^  a.  kc  -  ) 

“T  “  vT  ‘  Tu  '  Tf-i 

r.i  *  I.  Tl»  *  pt  '«  ,,-J  T.  |P«  ca  r  •  p|>  M>  Tb  ‘  Pc  >1  yj 


f  1 

(2)  —  <  a  <  Ta.  The  energy  balance  is 


^cond  1  ^cond  *  ^stored 


1  —  2  3  —  2 


jj) 

(^) 

(100) 


(101) 


f  >»' " 


(Taa-_a)  (Tm  ■  Ti)  +t"  (T’  •  T*)  =  [pa  ca  (ra  -  a) 


"1  Pb  cb  Tb  Pc 


Tc 
c  cc  2 


Tc|  (t2  -  T2) 


Solving  for  T2  , 


(102) 


Tm  k,  At  T;  kc  At  \-~)  r 

tc  -  -  lyr  1  t2  .  - 

Ta  P  .  '.i  fTa  aj  •  Pb  <-b  Tb  *  Pc  cc  y  tc  pn  t.,  ^Ta  -  nj  •  pb  cb  Tb  *  pt  cc  —  j  ' 


,  (Ta  -  a) 
kc  At  — 


‘I"-  ‘  a  (Ta  •  ■’)  *  •’!>  1  b  Tb  •  P.  •  c  y  I  Tc  {pa  S,  (Ta  *  »)  *  Pb  M,  Tb  ‘  Pc  «Y  y  j 


(103) 


w®-  ***■  - 


Heat  flow  in  the  second  thicn  iayer  is  determined  from  those  equa¬ 
tions  presented  in  Section  II  for  a  nonrecession  case. 


2.  Cylinder  Ablation -Conduction 

The  energy  balance  for  the  radial  heat  flow  toward  the  center- 
line  in  a  cylinder  is  basically  the  same  as  the  flat-plate  case,  with  the 
exception  that  average  areas  or  changing  areas  are  considered  for 
conduction  and  storage.  Once  the  ablation  front  reaches  ati  interior 
nodal  point,  the  ablation  distance  "a"  is  reduced  by  the  value  of 
Ta  (0  £  a  <  t a) ,  and  the  finite-difference  calculation  procc  ss  starts 
over  again.  In  addition,  the  radius  to  the  original  7;  fid:..'  point  is 
reduced  by  the  value  Ta. 


a.  General  Thick  with  6a  >  ra  (Figure  22) 
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or 
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=  Pa  ca  Ta 


At 


Az 
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whore  A  =  0L(R-/T2  + 


Ta  -  a 


)  -  2 


A  =  OL  [R  ->  t2 


T  a 


3-2 


A2  =  OL 


(«  -E  t.) 


Lotting 

Pa  a 


kdAt 


Pa  ca  Ta 

I 

and  solving  for  T2 


71  •  Z  =  R  -E  T* 


q  ,  A.  +  q  .  A  -  q  A2 

cond  i-2  cond  3-2  stored 
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1-2 
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A  =  0l(R£t 


3-  2 
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k  At 

Pa  ~  -r-2  ’  Z  =  R  'E  T*  * 

Pa  ca  Ta 

and  solving  for  T2i 


(108) 


Tm  Pj 


Tt 


>  — — r--\  . /  zTT^-ajfa-a) 


Pa 


IZ  +• 
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As  long  as  the  remaining  wall  thickness  of  the  material  undergoing 
recession  is  equal  to  one  or  more  thar  one  t,  the  interior  and  interface 
equations  presented  in  Section  II  are  used  to  solve  the  heat  balances 
throughout  the  structure  away  from  the  receding  surfa<  •  . 

b.  Special  Thick  with  6a  S  Ta  (Receding  Surface  <  Ta  from 
Backside)  (Figure  2i) 

Ta 

(1)  0  S  a  <  —  ,  The  energy  balance  for  T2  is 
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Pa  = 


ka  At 
Pa  ca  Ta2 


,  Z  =  R  "23  Tz  » 


and  solving  for  T2» 
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(2)  <  a  <  Ta.  The  energy  balance  for  T2  is 


^cond  j^2  “  %eti  Ai  “  ^stored 
1—2  2 


or 
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=  Pa  ca 
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Pa  ca  Ta 


and  solving  for  T 2 
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Lotting  7.  -  R  -  ^t2  and  solving  for  T 2  , 
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d.  Special  Thick- Thick  with  &a  <  ra  (Figure  25) 


(1)  0  <  a  <  —  .  The  energy  balance  for  T2  is 
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Letting  Z  =  R  -]£]t2  and  solving  for  T2, 
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(2)  —  <  a  <  Ta.  The  energy  balance  for  T2  is 


(124) 


^cond  j^2  +  c*cond  ~  ^stored  (A2  +  A,}) 


l—  2 


3—2 


(125) 


or 


A 

3-2 


T^-a)  (T™  -  *)  ,V  ^  (*  -  *) 

-jpa  ca  ^Ta  -  aj  A3  +  pb  cb -~  A4 


(T2  -  T2) 
At 


(126) 


50 


-Ml ■'»' 


where 
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Letting  Z  -  R  -]£  T2  and  solving  for  T2, 
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Heat  flow  in  the  second  thick  layer  is  determined  from  those 
equations  presented  in  Section  II  for  a  nonrecession  case. 

e.  Special  Thick-Thin-Thick  with  6a  £  Ta  (Figure  26) 
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Figure  26. 
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(i)  0  £  a  <  —  .  The  energy  balance  for  T2  is 

£ 
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I 
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(2)  —  <  a  <  Ta  .  The  energy  balance  for  T2  is 
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Heat  flow  in  the  second  thick  layer  is  determined  from  those  equa¬ 
tions  presented  in  Section  II  for  a  nonrecession  case. 


3.  Sphere  Ablation 'Conduction 

The  energy  balance  for  the  radial  heat  flow  toward  the  center 
in  a  sphere  is  basically  the  same  as  radial  heat  flow  in  a  cylinder  if 
the  average  areas  are  modified.  Appendix  A  of  Report  No.  RS-TR- 
65-13  gives  some  of  the  derivations  for  average  areas  for  spheres. 
Once  the  ablation  front  reaches  an  interior  nodal  point,  the  ablation 
distance  "a"  is  reduced  by  the  value  of  ra(0  <  a  <  ra),  and  the  finite- 
difference  calculation  process  starts  over  again.  In  addition,  the 
radius  to  the  original  Ti  nodal  point  is  reduced  by  the  value  ra. 


a.  General  Thick  with  6a  >  Ta  (Figure  27) 
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Figure  27. 


(1)  0  <  a<  —  .  The  energy  balance  for  Tz  is 
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(2)  —  <  a  <  Ta.  The  energy  balance  for  T2  is 
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As  long  as  the  remaining  wall  thickness  of  the  material  undergoing 
recession  is  equal  to  one  or  more  than  one  t,  the  interior  and  interface 
equations  presented  in  Section  II  are  used  to  solve  the  heat  balances 
throughout  the  structure  away  from  the  receding  surface. 


b.  Special  Thick  with  6a  <;  Ta  (Receding  Surface  s  Ta  from 


Backside)  (Figure  28) 
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(1)  Os  a<  —  .  The  energy  balance'  for  T2  is 
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(2)  —  <  a  s  Ta.  The  energy  balance  for  T2  i 
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or 
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c.  Special  Thick- Thin  with  6a  £  Ta  (Figure  29) 
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Figure  2Q. 


( 1 )  0  <  a  <  —  .  The  energy  balance  for  T2  is 
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(2)  ~  <  a  -  Ta*  The  energy  balance  for  T2  is 
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Letting  Z  =  R  -]£]t2  and  solving  for  T2, 
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(2)  —  <  a  <  t a •  The  energy  balance  for  Tz  is 
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Letting  Z  =  R  -Yl  tz  ,  and  solving  for  Tz  , 
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Heat  flow  in  the  second  thick  layer  is  determined  from  those  equa 
tions  presented  in  Section  II  for  a  nonrecession  case. 
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Section  IV.  HEAT  CONDUCTION  AFTER  ABLATION  TERMINATES  (T^ond  T*) 


With  the  equations  derived  for  heat  conduction  prior  to  ablation 
and  during  ablation,  the  heat  conduction  equations  after  ablation  stops 
must  be  considered.  The  equations  for  conduction  prior  to  ablation 
cannot  be  used  after  ablation  ceases  without  some  modification  to  the 
temperature  grid  because  there  is  no  assurance  that  the  ablation  into 
a  new  r  will  be  zero,  i.  e.  ,  the  ablation  usually  will  not  cease  with  the 
receded  surface  coinciding  with  an  original  temperature  node  location. 
Although  it  is  possible  to  select  a  new  temperature  grid  system  for  the 
material  left  after  surface  recession  ceases  and  to  obtain  proper  tem¬ 
perature  of  new  nodal  points  based  on  interpolations  from  the  calculated 
temperature  gradient  when  ablation  ceases,  an  approach  is  taken  herein 
whereby  the  original  grid  remains  unchanged.  Heat  conduction  equations 
for  T3  through  Tn  are  derived  in  the  same  manner  as  before  and  during 
ablation.  The  equations  used  to  calculate  Tj  and  T2  are  derived  in 
Paragraphs!,  2,  and  3. 


1.  Flat  Plat* 

a.  General  Thick  (Figure  16,  Tm  Taking  on  the  Value  Ti) 
with  6  a  >  t3 

Ta  i 

(1)  0  <  a<  —  .  The  energy  balance  for  li  is 


^netQ  ^  ~  ^cond  ^  “  ^stored  ^ 
1—2  1 


For  the  flat  plate  conduction,  A  may  be  taken  as  unity;  then, 


(164) 


Oneto  -  J—)  (Ti  -  Ti)  -  Pa  =a  (t  ‘  *)  AA  '  (165> 


The  energy  balance  for  Tz  may  be  taken  as 

^cond  ^  +  ^cond  A  =  ^stored  ^ 

1—2  3—2  2 


or 


(t;  '  t;) = p 


a  ca  Ta 


(Tz  -  T2) 


At 


(166) 


(167) 
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This  equation  is  the  same  as  Equation  (74)  while  the  surface  is  ablating, 
except  for  Tm  taking  on  the  value  T{. 


Since  T3  through  T^,  including  any  interfaces,  can  be  calculated 
by  using  ordinary  forward  finite-difference  methods,  Equations  (165) 
and  (167)  have  two  unknowns.  Solving  Equations  (165)  and  (167)  simul¬ 
taneously  and  letting 


^a  ^ _ 

Pa  ca  Ta2 


results  in 


,  l~^r~  *  p»  [■  ♦  ~V~]|  Ti  zrrz  1  r/S  vt  *  fT‘  * Ti  fo] 
,=  -#?•■)  lr,  \  |¥-*\ 

- ^ -  +  aa  1^1  +  QaJ  +  Pal. — 7^ — ) 


(Ta  '  >  l-r  -  a. 


Pa[  [T,  +  Pa  Tj]  +  Pa  ^TT  +  T,  Pa 

(Ta  -  a)  -  a)  p  1  It 

Ta  1 

(2)  —  <  a  <  t  .  The  energy  balance  for  Ti  is 
2  a 


(168) 


(169) 


9net0  ^  ~  ^cond  ^  "  ^stored  ^ * 

1  —  2  1 

For  the  flat- plate  conduction,  A  may  be  taken  as  umcy;  then, 


%et0  -  PThj  (t!  -  Tl)  ■ 


,neto  (ra-a)  V1  ' 

The  energy  balance  for  T->  may  be  taken  as 

^cond  ^  +  ^cond  ^  ~  ^stored  ^ 


(170) 


(171) 


(172) 


<T7h)  fa  • Tfa  £  fa  -  T>) 

-a  ^  -  a)  • 


(173) 
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This  equation  is  the  same  as  Equation  (78)  while  the  surface  is 
ablating,  except  for  Tm  taking  on  the  value  T{.  Solving  Equations  (171) 
and  (173)  simultaneously  and  letting 


= 


ka 

Pa  ca  Ta2 


gives 


%et. 


/Ja\ 
’  \ka/ 


(x  - a)  (T*  - a) 


*•  Pa  +  Pa 


(^) 


'3t, 


+  T*' 


+  tJ  oa 


3t„ 


+  Pa 


(174) 


and 


'  3t. 


-  a 


^neto  \  kaJ  +  Tz 


3t. 


-  a 


+  Pa 


+  T3  Pa 


(175) 


It  should  be  noted  that  the  denominators  of  Equations  (168)  and  (169) 
are  identical  as  are  the  denominators  of  Equations  (174)  and  (175). 

When  "a"  — ►  Ta,  Equations  (174)  and  (175)  for  T{  and  T2  reduce  as 
expected  to 


t;=t;  =q„et0^  +  Ti 


(176) 


b.  Special  Thick  (Receding  Surface  <  t3  from  Backside, 
Figure  18,  Tm  Taking  on  Value  Tj)  with  6a  <  Ta 


When  the  receding  surface  is  less  than  one  incremental 
t  from  another  material  surface  or  backside,  special  considerations 
must  be  made  for  all  material  combinations  normally  experienced. 

Ta 

(1)  0  <  a<  —  .  The  energy  balance  for  Tj  is 


^neto  ^  ~  ^cond 


^  ~  ^stored  ^ 


(177) 


For  the  flat-plate  conduction,  A  may  be  taken  as  unity;  then, 


qnet0  ‘  |T&  _  a)  (Tl  '  T?')  =  Pa  Ca  (  2  "  a) 
The  energy  balance  for  T2  may  be  taken  as 
qcond  ^  '  qnet^  ^  ~  ^stored  ^ 


ir'i  -  Ti) 
At 


(178) 


(179) 


(Ta  -  a)  (Tl  ‘  Tz)  '  qneti  =  Pa  °a  2 


(Tz  -  T2) 
At 


(180) 


Solving  Equations  (178)  and  (180)  simultaneously  and  letting 
ka  ^ 

Pa  q  „  2 

Pa  ca  Ta 


gives 


Nr  +  N  [T>  NN  *  5rfrr  +  ^  [T*  -  4^ 

(t  * a)  (Ta " a)  .  '  It  -  A 


+  Pa  +  2fla 


(181) 


/Ta  a\  1 

/T  a  \/  \ 

1 

T  (_2 '  .  qm!t°  At 

x 

(t  ■  #a  ■  a) ,  n 

2  qnetj  ^ 

.  \  Ta  /  Pa  Ca  Tn  , 

T 

L  Ta*  +  Pa 

l  2  Pa  Ca  Ta  j 

/Ta  \/  \  /  Tn  \ 


" - 1  +Pa+2Pa  -7- 


(182) 


*  3  i 

(2)  —  <  a<  Ta.  The  energy  balance  for  Tj  is 

9  <  A-q  ,  A  =  q  .  ,  A 

lnet0  'cond  'stored 

1  —  2  1 

For  the  flat-plate  conduction,  A  may  be  taken  as  unity;  then 


'  ''  '^' *  './‘v  k<fy  ^ /;!;%,  'V\  *'\  "  -  ■  *~ 


'<*’  i 


The  energy  balance  for  T2  is 


^cond  A  "  %et;  A  =  ^stored  A 
1— *2  , 


(185) 


or 


|^y  (« -  ti) p.  e a(I^)^.  < 

Solving  Equations  (184)  and  (186)  simultaneously  and  letting 
ka  At 


Pa  ca  Ta2 


results  in 


t;  =  t2  + 


Ta 


(Ta  -  a) 

^netj  At 


pa  ca  Ta 


qnet°(k!)  Pa+(~ 


(187) 


and 


T'z  =  T2  4 


4.  Ta  f„  .  /Ii\  „  8 net,-  At 

'  (r.  -  a)  [<■«*„  (kJ  Pa  -  7r^Ta 


(188) 


C-  Special  Thick-Thin  (Figure  19,  T„  Taking  on  Value  T*) 

unf  h  A  1  •  1/ 


witk  5a  <  Ta 


(1)  0  <  a  <  ~  .  The 


energy  balance  for  T[  is 


qnetQ  A  '  ^cond  A  =  ^stored  A 


1  —  2 


(189) 


For  the  flat- plate  conduction,  A  may  be  taken  as  unity;  then, 

q-0  -  prh)  (Ti  -  Ti)  1  ^  ca(r  ■  a)  •  (190) 
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The  energy  balance  for  T'z  may  be  taken  as 
9cond  A  -  ^neti  ^  =  ^stored  ^ 


1  —  2 


(191) 


or 


(r “Vl)  (Tl  -  Tz)  ‘  qneti  =  (Pa  ca  2- 

\  (Tz  -  T2)  . 

+  Pb  cb  Tb)  ~st — 


(192) 


Solving  Equations  (190)  and  (192)  simultaneously  and  letting 

Pa  = 


ka  At 


Pa  ca  Ta 


gives 


T! c 


k.  At 


*  7  ra  \ 

Ta(Pa  Ca"T  *  Pb  <-b  Tbj 


ri(i_:V  Pl 

1  \  Ta  J  Oa  ca  Ta 


T,  -- 


Pncti  At 


Pa  ca  -j*  ♦  Pb  cb  Tb 


(t4  •  »)(t  '  a) 


ka  At’ 


♦  P;, 


SI 

\  Ta 


.  _ 

Ta  Pa  Ca  —  '  Pb  'b  >b 


(193) 


and 


ka  At 


>  (pa  ca  Y  1  fb  cb  Tb) 


.1  [t/xh)  *  assail 
)]  v  TJ  •  c* Ta 


(7»  •  4t  • *) 


*  Pa 


T,  -- 


9netj  ' 


(Tj  ~  j)(t  '  a) 


'  Pa  * 


ra(pa  ca-^  ♦  Pb  'b  Tb) 


(194) 


(2) 


y  <  a  <  Ta. 


The  energy  balance  for  Tj  is 


qnetQ  ^  “  qcond  ^  ~  ^stored  ^ 


1—  2 


For  the  flat- plate  conduction,  A  may  be  taken  as  unity;  then, 


w„  -  fzh)  (Tl  - Ti)  ■ 0 


(195) 


(196) 
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The  energy  balance  for  T  is 


9cond  A-  ~  ^neti  ^  =  ^stored  ^ 


(197) 


f^rh)  (t!  •  Tl)  -  q„etl=[paCa(Ta-  a) 


+  Pb  c  b  T  b 


(t'z  -  TZ) 


(198) 


Solving  Equations  (196)  and  (198)  simultaneously  gives 


“  ^2  +  ^lnet. 


|Pa  ca  (Ta  -  a)  +  Pb  cb  Tb|  k 


Ta  -  a| 


9netj  At _ 

jPa  ca  (ra  -  a)  +  pb  cb  xb| 


(199) 


At (9net0  ‘  9netJ 


T,  =  T,  +  - 1 - - j 

ra  ca  (Ta  -  a)  +  Pb  cb  Tb} 


(200) 


d.  Special  Thick-Thick  (Figure  20,  Tm  Taking  on  Value  T{  ) 
with  6a  <  t a 

(1)  0  5  a  <  y  .  The  energy  balance  for  Tj  is 


9net0  ^  ~  ^cond  ^  ~  9stored  ^ 
1  —  2  l 


(201) 


ka  (TX  -  T z) 
(Ta  "  a)  ' 


=  Pa  ca 


-  Ti) 
At 


(202) 
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t  , 

> 


The  energy  balance  for  T2  is 

^cond  ^  +  ^cond  A  =  ^stored  ^ 
1—2  3—2  2 


jrhf  H*  £  (Ti-Ti)'hc-T 

.  .  .  _  .  .  _  Tcl  (t2  -  Tz) 


+  Pb  cb  Tb  +  Pc  cc  ~2 


Solving  Equations  (214)  and  (216)  simultaneously  and  letting 


B  -  Pa  ca  T  +  Pb  <=b  Tb  +  Pc  cc  ~y~ 


ka  At 

[3a  =  - - - -2 

Pa  ca  Ta 


ka  At 

P»  *  av. 

Ta  +  t  c 

l  Ta  ) 

LO  t/jluV  +  (ia [t2  +  TJ 44rl 


(2)  —  <  a  <  ra.  The  energy  balance  for  Tj  is 


or 


^notQ  ^  ~  ^cond  ^  =  ^stored  ^ 
1  —  2  1 


q""°  -  prr^  (Ti  ■  Ti)  ■ 0 

The  energy  balance  for  T2  is 
^cond  ^  +  ^cond  ^  =  ^stored  ^ 

1—2  3—2  2 


(219) 


(220) 


(221) 


or 


fra  -  a)  (T^  ‘  T^)  +  1“  (T*  ‘  T*)  =[pa  ca  (Ta  -  aj  +  pb  cb  rb 

:  J  (t2  ~  TZ)  . 


*  Pc  cc  2 


(222) 


Solving  Equations  (220)  and  (222)  simultaneously  and  letting 


B  =  Pa  ca 
results  in 


(t3  -  a) 


+  Pb  cb  Tb  +  pc  cc 


ka  At  kc  At  /ra  -  a\ 

qnc,o  (ka) 

Ta  *  a  .  Ta  *c  1  'a  / 

T  kc  At 

.  Ta  D 

*  Tl  +  Ts  7Tb" 

(223) 


i  ♦ 


kc  At 


and 


Ti  =■ 


^neto  ^  rr,  m'  ^C  At 

- 2 -  +  Tz  +  T3  - — 

_ Tq  D 


B 


1  + 


kc  At 

TC  B 


(224) 
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The  energy  balance  for  T2  is  the  same  as  Equal 
replaced- by  T{. 


1.0.8)  if  Tm  is 


Solving  Equation  (230)  and  modified  Equation  (108)  simultaneously 
and  letting 


Pa 


ka  At 
Pa  ca  Ta2 


and  Z  =  R  -£t2  results  in 


and 


Tz 


(232) 


where 


Ta  ”  a 

Bj  =  Z  +  - 

Ta 

B3  =  z  -  T 

— 

B4  =  z+i_ 


i- 


1 


Bs  -  R  -  a  . 


b.  Special  Thick  (Exposed  Surface  <  t  from  Backside) 
(Figure  23,  Tm  Taking  on  Value  Ti)  with  6a  .■S'  ra 

Ta  t 

(1)  0  £  a  <  —  -  .  The  energy  balance  for  Tj  is  the  same 

as  Equation  (226). 

The  energy  balance  for  T2  is  the  same  as  Equation  (111)  if  Tm  is 
replaced  by  T(. 

Solving  Equation  (226)  and  modified  Equation  (111)  simultaneously 
and  letting 

ka 

3  = - 

Pa  ca  Ta2 

and  Z  =  R  -  ^]t2  gives 


where 


„  „  Ta  - 

B*  =  Z  + - 2 

B2  Z 

ila 

Z  < 


B5  R  -  a 

B6  Z  <  y 
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a  1 

(2)  —  <  a  <  Ta.  The  energy  balance  for  Ti  is  the  same 

as  Equation  (230). 

The  energy  balance  for  T2  is  the  same  as  Equation  (114)  if  Tm  is 
replaced  by  T{. 

Solving  Equation  (230)  and  modified  Equation  (114)  simultaneously 
and  letting 

ka  At 

^  ^  O  C  T  2 

pa  '-a  a 


and  Z  =  R  -^t2  results  in 


Ti  =  T2  +  7 


F  - a) 

—  I^A 

7  a  \Bi/ 


qnetj  At  /B2 
pa  ca 


°Inet0 
1 


Ua /  iBi/  {  a  \  Ta  /  j 


(235) 


and 


T2  =  T2  + 


(Ta-a)  |_ 

-  qneti _  /BA 

a  '-a  Ta  \Bi/ 


^netQ  Pa 


0© 


Pa  C 


(236) 


where 

Ta  ■  a 

B1  =  Z+— — 

B2  =  z 

B5  =  R  -  a  . 

c.  Special  Thick-Thin  (Figure  24,  Tm  Taking  on  Value  Ti) 
with  6a^  t a 

Ta  1 

(1)  0  £  a  <  — .  The  energy  balance  for  Ti  is  the  same  as 

Equation  (226). 
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The  energy  balance  for  T2  is  the  same  as  Equation  (117)  if  Tm  is 
replaced  by  T(. 


Solving  Equation  (226)  and  modified  Equation  (117)  simultaneously 
and  letting 

ka  At 


and  Z  =  R  2  gives 


whe  re 

®  +  ”4”)  ca  ~2~  +  Pb  cb  Tb 


Bi  =  Z 


Ta  *  a 
2 

^Ta 


-  a 


B4  =  Z  + 


B5  -  R  -  a 
B7  -  Z  -  Tb  . 


»  ^  I 

(2)  —  <  a  <  t3.  The  energy  balance  for  Tj  is  the  same 
as  Equation  (230). 

The  energy  balance  for  Ti  is  the  same  as  Equation  (120)  if  Tm  is 
replaced  by  Ti. 
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Solving  Equation  (230)  and  modified  Equation  (120)  simultaneously 
and  letting  Z  -  R  -  Y'.  T2  results  in 


Ti  r  Tz  +  %et. 


-f Hxl 


.Ml2)  (239) 


'  qnoti  m  l  B 


and 


.  Al  ^neto  (B5)  ‘  ^neti 

t2  =  t2  + - - - — g - 1 - 


(240) 


where 

B 


--  (Z  4  --~2— -)  Pa  ca  (Ta  ‘  a)  +  (Z  ’  t)  °b  Cb  Tb 


Bi  =  Z  + 


Ta  '  a 


B5  -  R  -  a 
B7  -  Z  -  rb 

d.  Special  Thick-Thick  (Figure  25,  Tm  Taking  on  Value  Ti) 
With  6a<  t a 

T3  through  Tn  are  calculated  from  the  general  heat  conduc¬ 
tion  equations  for  cylinders. 


(1)  0  <  a  <  ~~  .  The  energy  balance  for  Tj  is  the  same 

as  Equation  (226). 

The  energy  balance  for  T2  is  the  same  as  Equation  (123)  if  Tni  is 
replaced  by  Tj. 


Solving  Equation  (226)  and  modified  Equation  (123)  simultaneously 
and  letting  Z  R  -  7  .  t2  results  in 


(241) 


S3 


and 


t1 


(242) 


where 


Pa  ca 


+  Z 


Tb\  Tb 

T/Pb  CbT 


B!  =  Z  + 


Ta  "  a 


3t. 


-  a 


B4  -  Z  + 


B5  -  R  -  a 

Tb 

B«  =  z  -  T  • 

(2)  <  a  <  ra.  The  energy  balance  for  Ti  is  the  same 

as  Equation  (230). 

The  energy  balance  for  T2  is  the  same  as  Equation  (126)  if  Tm  is 
replaced  by  Tj. 


Solving  Equation  (230)  and  modified  Equation  (126)  simultaneously 
and  letting  Z  =  R  -£}t2  gives 


a  MM 

r.i  -  A 

k  .  At  kb  /~A  -  AV 

,,  <».*'  tb  <».l  (  ra  ) 

kb  At  <B*) 

(B) 

T’  (B) 

u>.) 

'  '  *.,(»> 


(243) 


and 


T2 


+ 


T2  +  T'3 


kb  At  (B8) 
"  b  (B) 


1 


kb  At  (Bg) 
+  _Tb  (B) 


(244) 
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whore 


B5  -  R  -  a 

„  „  Tb 

B8  -  Z  -  T  • 


c.  Special  Thick- Thin- Thick  (Figure  26,  Tm  Taking  on 

Value  Tj)  w  ith  6 a  <  Ta  '  ~ 

T3  through  Tj,  are  first  calculated  using  the  forward 
finite-difference  equations  for  cylinders. 

T  a  , 

(1)  0  <  a<  —  .  The  energy  balance  for  Tj  is  the  same 

as  Equation  (226). 

The  energy  balance  for  T2  is  the  same  as  Equation  (129)  if  Tm  is 
replaced  by  T{. 


Solving  Equation  (226)  and  modified  Equation  (129)  simultaneously 
and  letting 

k  At 


and  Z  =  R  t2  results  in 


) 

t 


J*- 

f. 


where 


B 


Ta 

Pa  ca  2 


Pb  cb  Tb  + 


(z 


Tc 

Pc  cc  2 


Ta  "  a 
Bi  =  Z+— — 

B4  =  Z  +  -~- 


B5  =  R  -  a 


Ta  1 

(2)  —  <  a  <  Va.  The  energy  balance  for  Ti  is  the  same 

as  Equation  (230). 

The  energy  balance  for  T2  is  the  same  as  Equation  (132)  if  Tm  is 
replaced  by  Ti. 

Solving  Equation  (230)  and  modified  Equation  (132)  simultaneously 
and  letting  Z  -  R  -]>^T2  gives 


and 


/B5\  ^  ,  kc 

qnot0  At  (— J  +  Tz  +  T3  — 


1  + 


kc  At 


whe  re 


B 


Pa  c*  lTa 


(T- - a) +  (z  -  t) 


Pb  Cl'  T 


b  Tb 


+  Z 


tc 

Pc  cc  2 


(248) 
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Ta  -  a 
Bl  =  Z+~r- 

B5  =  R  -  a 

Tc 

Bg  -  Z,  -  Tjj  -  . 


3.  Sphere 


a.  General  Thick  (Figure  27,  Tm  Taking  on  Value  Ti)  with 


6a  >  Ta 


T3  through  T|j,  including  any  interfaces,  are  calculated 
from  the  general  heat  conduction  equations  for  spheres. 

Ta 

(1)  0  <  a<  —  .  The  energy  balance  for  T{  is 


^neto  Al  “  ^cond  Az  "  ^stored 
1-*-  2  1 


or 


^netf 


Al '  fTT)  (Tl '  Ti)  Al 

/la  \(t!  -  Ti) 
V  2  /  At 


=  A3  pa  c 


where 


^  =  <}.  [(R  -  a)2] 


A2  =  4. 


A3  =  <j> 


3Ta 


(Ta  -  ») 


12 


£ 


-  a 


R  - 7  T  2  + 


JH 

12 


(249) 


(250) 


The  energy  balance  for  T2  is  the  same  as  Equation  (135)  if  Tm  is 
replaced  by  Ti. 
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Solving  Equation  (250)  and  modified  Equation  (135)  simultaneously 
and  letting  Z  r  R  -]T]t2  results  in  Equations  (227)  ana  (228)  where  the 
following  parameters  take  on  the  new  values  of 


(251) 

(252) 

(253) 

(254) 

(255) 


(2) 


a  <  Ta. 


The  energy  balance  for  is 


^net0  ^  “  ^cond  -^2  ~  ^stored 
1—2  l 


(256) 


or 


qnet0  A'  +  (Ti  '  Ti)  Az  =  0 


(257) 


where 


Ai  =  4>  (R  -  a) 
A2  =  <j> 


(Ta  -  *)' 


12 


The  energy  balance  for  T2  is  the  same  as  Equation  (138)  if  Tm  is 
replaced  by  Tj. 


Solving  Equation  (257)  and  modified  Equation  (138)  simultaneously 
and  letting  Z  =  R  -]T^t2  gives  Equations  (231)  and  (232)  where  the 
following  parameters  take  on  the  new  values  of 


88 


B5  =  (R  -  a)2  . 


(258) 

(259) 


(260) 

(261) 


b.  Special  Thick  (Exposed  Surface  <  Ta  from  Backside) 
(Figure  28,  Tm  Taking  on  Value  Tj)  with  6a  <  Ta~ 


•  3  I 

(1)  0  <  a<  —  .  The  energy  balance  for  Ti  is  the  same 

as  Equation  (250). 

The  energy  balance  for  T2  is  the  same  as  Equation  (141)  if  Tm  is 
replaced  by  T|. 


Solving  Equation  (250)  and  modified  Equation  (141)  simultaneously 
and  letting 

ka  At 

^  a  a  a  t  ^ 

Pa  ca  Ta 


and  Z  =  R  _23T2  results  in  Equations  (233)  and  (234)  where  the  following 
parameters  take  on  the  new  values  of 


(262) 

(263) 

(264) 

(265) 

(266) 


*4  7 


(2)  — -  <  a  <  ra.  The  energy  balance  for  Ti  is  the  same 
as  Equation  (257). 

The  energy  balance  for  T2  is  the  same  as  Equation  (144)  if  Tm  is 
replaced  by  Tj. 


Solving  Equation  (257)  and  modified  Equation  (144)  simultaneously 
and  letting 


Pa 


ka  At 
Pa  ca  Ta2 


and  Z  =  R  tz  results  in  Equations  (235)  and  (236)  where  the  following 
parameters  take  on  the  new  values  of 


(267) 

(268) 


Bs  =  (R  -  a)2  . 


(269) 


c.  Special  Thick-Thin  (Figure  29,  Tm  Taking  on  Value  T{) 
with  6a  <  Ta 


1  S  1 

(1)  0  <  a<  —  .  The  energy  balance  for  Tj  is  the  same 

as  Equation  (250). 

The  energy  balance  for  T2  is  the  same  as  Equation  (147)  if  Tm  is 
replaced  by  T{. 


Solving  Equation  (250)  and  modified  Equation  (147)  simultaneously 
and  letting  Z  =  R  -T!  T2  gives  Equations  (237)  and  (238)  where  the 
following  parameters  take  on  the  new  values  of 


Tbzj 

+  l2j  pbCbTb 


(270) 
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Bs  =  (R  -  a)2 
B7  =  (Z  -  rb)2 


(271) 


(272) 

(273) 

(274) 


cl  I 

(2)  —  <  a  <  Ta.  The  energy  balance  for  Ti  is  the  same 


as  Equation  (257). 


The  energy  balance  for  T2  is  the  same  as  Equation  (150)  if  Tm  is 
replaced  by  Ti. 


Solving  Equation  (257)  and  modified  Equation  (150)  simultaneously 
and  letting  Z  =  R  -'jT,  t2  gives  Equations  (239)  and  (240)  where  the  following 
parameters  take  on  the  new  values  of 


(275) 

(276) 

(277) 

(278) 


d.  Special  Thick- Thick  (Figure  30,  Tm  Taking  on  Value  Tj) 
with  6a  <  t a 


T3  through  Taare  calculated  from  the  general  heat  conduc¬ 
tion  equations  for  spheres. 


(1) 

Equation  (250). 


0  <  a  < 


The  energy  balance  for  Ti  is  the  same  as 
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Tin*  energy  balance  for  T2  is  the  same  as  Equation  (153)  if  Tm  is 
replaced  by  Ti. 


Solving  Equation  (250)  and  modified  Equation  (153)  simultaneously 
and  letting 

k.  At 

S-  = - - - — 

Pa  ca  Ta 2 


and  Z  =  R  -]T}t2  results  in  Equations  (241)  and  (242)  where  the  following 
parameters  take  on  the  new  values  of 


(279) 

(280) 

(281) 

(282) 

(283) 


(2) 


a  <  T' 


as  Equation  (257). 


The  energy  balance  for  Tj  is  the  same 


The  energy  balance  for  T2  is  the  same  as  Equation  (156)  if  Tm  is 
replaced  by  Tj. 


Solving  Equation  (257)  and  modified  Equation  (156)  simultaneously 
and  letting  Z  --  R  t2  gives  Equations  (243)  and  (244)  where  the  following 
parameters  take  on  the  new  values  of 


'12 


(284) 


(285) 


Bs  =  (R  -  a)2 


(286) 


B8  = 


(287) 


e.  Special  Thick-Thin-Thick  (Figure  31,  Tm  Taking  on 
Value  Tj)  with  6a  £  ra 


T3  through  T„  are  first  calculated  using  the  forward 
finite-difference  equations  for  spheres. 

Ta  1 

(1)  0  <  a  <  —  .  The  energy  balance  for  Ti  is  the  same 

as  Equation  (250). 

The  energy  balance  for  Tz  is  the  same  as  Equation  (159)  if  Tm  is 
replaced  by  T{. 

Solving  Equation  (250)  and  modified  Equation  (159)  simultaneously 
and  letting 

ka  At 

6-  = - ■, 

Pa  ca  Ta 


and  Z  =  R  2  results  in  Equations  (245)  and  (246)  where  the  following 
parameters  take  on  the  new  values  of 
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Bs  =  (R 

[/ 

-  a)2 

M2  Tc2| 

(291) 

b9  =  [p 

-Tb-T/  +T2  J  * 

(292) 

(2)  -£■  <  a  <  Ta.  The  energy  balance  for  T{  is  the  same 
as  Equation  (257). 

The  energy  balance  for  T2  is  the  same  as  Equation  (162)  if  Tm  is 
replaced  by  T{. 


Solving  Equation  (257)  and  modified  Equation  (162)  simultaneously 
and  letting  Z  =  R  "J]JT2  resuRs  Equations  (247)  and  (248)  where  the 
following  parameters  take  on  the  new  values  of 


(293) 

(294) 

(295) 

(296) 
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Section  V.  CRITERIA  TO  STOP  ABLATION  OR  SURFACE  RECESSION 


Once  a  material  surface  has  reached  the  critical  temperature  for 
ablating,  melting,  or  subliming,  and  the  surface  starts  to  recede,  the 
criteria  for  stopping  the  ablation  (recession)  must  be  established. 

Under  normal  conditions  the  heating  rate  to  the  ablating  surface  is 
decreasing  with  time  when  the  ablation  stops.  With  this  in  mind,  it 
was  decided  to  examine  the  Ti  equations  applicable  to  post- ablation 
heat  flow  to  see  if  a  criterion  could  be  established.  Since  qneto  the 
only  driving  parameter  found  in  the  Tj  equations  in  Section  IV,  there 
are  critical  values  of  qnetD  bel  n  which  the  exposed  surface  tempera¬ 
ture  cannot  be  maintained  at  the  ablating  temperature.  That  is,  more 
heat  is  being  conducted  internally  from  the  heated  surface  than  is 
available  at  the  heated  surface.  The  following  list  shows  how  to  find 
the  net  heating  rate  at  which  ablation  ceases  for  all  structural  arrange¬ 
ments  considered  in  this  report.  In  each  equation  listed,  Ti  is  set 
equal  to  Tni  before  solving  for  the  critical  value  of  qnet0- 


1.  Flat  Plate 

a.  General  Thick 

(1)  0  <  a  <  ^  . 

Equation  (168)  solved  for  qnetQ*  (297) 

(2)  <  a  <  Ta. 

Equation  (174)  solved  for  qaet0*  (298) 

b.  Special  Thick 

(1)  0  *  a<  y  . 

Equation  (181)  solved  for  qnetQ  .  (299) 

Ta 

(2)  ~  <  a  <  Ta. 

Equation  (187)  solved  for  qnet0'  (300) 
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Special  Thick-Thin 


(1)  0  <  a  <  ~  . 

Equation  (193)  solved  for  qncto. 

(2)  —  <  a  <  Ta. 

Equation  (199)  solved  for  qnct0> 

d.  Special  Thick- Thick 

(1)  0<a^  . 

Equation  (205)  solved  for  qnetD- 
(2>  ~  <  a  ^  Ta- 
Equation  (211)  solved  for  qnet  . 

e.  Special  Thick- Thin- Thick 

(1)  0  <  a  <  ~  . 

Equation  (217)  solved  for  qnet0* 

(2)  —  <  a  £  Ta. 

Equation  (223)  solved  for  9nct0- 

2.  Cylinder 

a.  General  Thick 

(i)  0  <  a  <  . 

Equation  (227)  solved  for  qneto- 


(301) 


(302) 


(303) 


(304) 


(305) 


(306) 


(307) 
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(308) 


(2)  <  a  ^  Ta- 

Equation  (231)  solved  for  qnet0  • 

b.  Special  Thick 

(1)  0<  a<~  . 

Equation  (233)  solved  for  qnet0* 
Ta 

(2)  -j-  <  a  <  Ta. 

Equation  (235)  solved  for  qnet0* 

c.  Special  Thick- Thin 

(1)  Q-  a-  T  • 

Equation  (237)  solved  for  qnet0’ 

(2)  —  <  a  <  Ta. 

Equation  (239)  solved  for  qnet0  • 
\  d.  Special  Thick-Thick 

(1)  0  <  a<  -£■  . 

Equation  (241)  solved  for  qm»t  • 
Ta 

(2)  —  <  a  <  Ta. 
Equation  (243)  solved  for  qnetQ- 


J 

(309) 


(310) 


(311) 


(312) 


(313) 


(314) 
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e. 


Special  Thick- Thin-Thick 


(1)  0  <  a  <  —  . 


Equation  (245)  solved  for  qnet0*  (315) 

(2)  -j-  <  a  <  Ta. 

Equation  (247)  solved  for  qnet0*  (316) 


3.  Sphere 

a.  General  Thick 

(1)  0  <  a<  . 

Solve  for  qnetQ  in  Equation  (227)  with  Equations  (251) 

through  (255)  included.  (317) 

(2)  --<*<,  Ta- 

Solve  for  qnetQ  in  Equation  (231)  with  Equations  (258) 

through  (261)  included.  (318) 

b.  Special  Thick 

(1)  0  <  a  <  —  . 

Solve  for  qnct0  in  Equation  (2  33)  with  Equations  (262) 
through  (266)  included.  (319) 

'  (2)  -£  <  a  <  ts. 


Solve  for  qnet0  in  Equation  (235)  with  Equations  (267) 
through  (269)  included. 
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(320) 


c.  Special  Thick-Thin 

(1)  0  ^  a  -  T  ‘ 

Solve  for  qnet0  in  Equation  (237)  with  Equations  (270) 

through  (274)  included.  (321) 

(2)  y  <  a  <  ra. 

Solve  for  qnet0  in  Equation  (239)  with  Equations  (275) 

through  (278)  included.  (322) 

d.  Special  Thick- Thick 
(D  0  £  a<;  ~  • 

Solve  for  9netQ  in  Equation  (241)  with  Equations  (279) 

through  (283)  included.  (323) 

(2)  -^  <  a  <  Ta. 

Solve  for  qnet0  in  Equation  (243)  with  Equations  (284) 

through  (287)  included.  (324) 

e.  Special  Thick-Thin-Thick 

(1)  0  <  a  <y  • 

Solve  for  qnetQ  in  Equation  (245)  with  Equations  (288) 

through  (292)  included.  (325) 

(2)  —  <  a  <  Ta. 

Solve  for  qncto  in  Equation  (247)  with  Equations  (293) 

through  (296)  included.  (326) 
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Section  VI.  CONCLUSIONS 


Excellent  agreement  is  obtained  between  data  generated  by  the 
combined  forward- backwa rd  finite-difference  equations  and  selected 
exact  analytical  equations  for  a  flat  plate  undergoing  surface  recession, 
when  proper  time  and  distance  increments  are  chosen  in  relation  to 
material  thermal  properties  and  surface  recession  rates.  The  com¬ 
bined  forward- backward  finite-difference  ablation- conduction  method 
is  most  accurate  when  the  amount  of  material  removed  in  a  calculation 
time  increment  is  equal  to  or  less  than  one-fourth  of  the  selected 
incremental  distance  between  temperature  nodes. 

Using  the  methods  presented  in  this  report,  simultaneous  conduc¬ 
tion  and  ablation  calculations  for  transient,  radial  heat  flow  in  spheres 
and  cylinders  and  one- dimensional  heat  flow  in  flat  plates  requires  a 
negligible  increase  in  computer  time  over  a  nonreceding  case  with  all 
other  parameters  being  identical.  This  is  true  primarily  because  the 
majority  of  the  equations  used  are  identical  with  those  used  in  nonreces¬ 
sion  cases. 

Forward-backward  finite-difference  equations  can  be  mixed  to 
achieve  simplicity  and  to  avoid  instability  in  equations  for  heat  flow 
near  the  surface  of  ablating,  subliming,  or  melting  structural  materials. 

Centripetal  ablation  and  heat  flow  equations  for  cylind 
spheres  can  be  modified  by  minor  sign  changes  to  obtain  ce. 
ablation  and  heat  flow  equations  for  cylinders  and  spheres. 
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Appendix  A 


DERIVATION  OF  CENTRIFUGAL  HEAT  CONDUCTION  AND  ABLATION 
EQUATIONS  AND  A  COMPARISON  OF  THESE  EQUATIONS  WITH 
THOSE  FOR  CENTRIPETAL  HEAT  FLOW 

1.  Cylinder 

«s 

The  energy  balance  for  the  radial  heat  flow  away  from  the 
centerline  (centrifugal)  of  a  cylinder  (Figure  22)  will  now  be  derived 
for  the  General  Thick  case  during  ablation.  For  centrifugal  flow,  R 
is  the  inside  radius  instead  of  the  outside  radius  used  for  centripetal 
flow  in  Sections  I,  III,  IV,  and  V  of  this  report. 

Ta 

a.  0  <  a  <  — 

The  energy  balance  for  T2  is 

qcond  j^2  +  ^cond3^2  "  ^stored  ^2  (327) 

i-~  2  3—2  2 


-v  L-  Ti)  A  +  ~  (T3  -  Ti)  A 

(Ta  -  a)  V  /  i-z  ' a  V  I  3-2 

(t 2  -  Tz) 


pa  ca  Ta 


where 


A  OL  R  +)  t ? 

1-2  LmJ 


A  -  OL  R  4 

3-2 


Az  -  OL 


(r  +£t2]  • 


Letting 


ka  At 


a  p  C  T  2 

pa  a  a 


Z  =  R  +  £' 


(328) 

(329) 

(330) 

(331) 
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and  solving  for  T2 


Ti 


T„,  Pa 


Ta  *  a 

7  *  Ta 

?.  -  2  - 

♦  Ti  Pa 

X  *  2 

,  T  /"a*  a) 

_ _ s _ _ 

Z 

l  Ta  / 

Ti.L  iaj 

Pa 

b. 


Ta  "  a 

i  Ta 

2 

•  Pa 

Y‘ '  i 

/Ta  -  a\ 

•Z  , 

L  z 

1  Ta  ) 

Ta 

2  <  a 

<  Ta 

T  a 


The  energy  balance  for  T2  is 
^cond  j^2  +  ^cond  3^2  ~  ^stored 


1  —  2 


3  —-2 


or 


(T™ '  ri)  A  +  $  (Ti  - T5)  A 

(t;  -  t.) 


pa  ca  l  2 


At 


Az 


where 


A  =  GL 
1-2 


A  =  GL 
3-2 


A2  -  GL 
Letting 


♦E 


R  +>  r2  - 


Ta  -  a 


R  + 


V'  Ta 

2Jz+t 

"lb 


Ta  -  a 


®a  = 


ka  At 


0  C  T 

*a  a  a 


2  * 


Z  =  R  +  £t2 


and  solving  for  T2 


i.  >, 


M  IrM  • " &  #' J 


'  V 

7  •  y 

(V  •  •)( 

7  ^ 

4V  J 

7  ^ 

(332) 


(333) 


(334) 


(335) 

(336) 

(337) 


(  3  38) 
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Equations  (106)  and  (332)  are  similar  as  are  Equations  (109)  and 
(338).  Equation  (332)  can  be  obtained  from  Equation  (106)  by  substituting 
a  multiplication  factor  into  Equation  ( 106).  The  same  is  true  for  Equation 
(109)  to  obtain  Equation  (338).  The  resulting  substitutions  and  equations 
are: 


Z  =  R  -  (XQ)£t2 


(339) 


c.  0  <  a  <  -j- 


r  -  a  1 

T 

7  *  (XQ)  —~y~ 

/  -  <XO)T 

7.  J 

7.  J 

\  J  Tz\ 

t.  -  a 

7.  4  (XQ) 

7.  -  <XO)~ 

/Ta  ‘  *\ 

(  Ta  -  a 

7. 

'  Pa 

7. 

\  TJ 

Ta 

(340) 


d.  — -  <  a  <  Ta 


12 


/.  *  (XQ) ---j— 


7  ♦  (XQ)- 


2a 


7  ■  (XQ>~ 


7  •  (XQ)  - 


(xjJM 


f,  ■  a 
Z  •  (XQ)-~-- 


7  >  <XQ)- 


t,  •  2a 


*  H 


ft) 


7.  -  (XQ)  -f 


7  .  (XQ)- 


(x  ■*)(’*■ a) 


(341) 

Factor  XQ  takes  on  the  value  of  +1  for  centripetal  flow  and  -1  for 
centrifugal  flow.  The  appropriate  radius  is  that  radius  to  the  original 
Ti  nodal  point.  Only  one  set  of  equations  is  compared  here;  however, 
all  of  the  equations  for  the  cylindrical  flow  have  been  investigated  and 
can  be  arranged  in  this  form. 


Due  to  the  nondimensional  terms  used,  it  is  also  possible  to  use  a 
-R  value  for  R  and  achieve  the  same  results.  In  other  words,  the 
centripetal  equations  presented  in  the  main  body  of  this  report  can  be 
used  in  their  present  form  to  determine  centrifugal  heat  flow  results 
simply  by  inserting  the  radius  to  the  initial  inner  surface  as  a  negative 
value. 


2.  Sphere 

The  energy  balance  for  the  radial  flow  away  from  the  center 
in  a  sphere  (Figure  27)  will  now  be  derived  for  the  General  Thick  case 
during  ablation.  For  centrifugal  flow,  R  is  the  inside  radius  instead 
of  the  outside  radius  used  for  centripetal  flow  in  Sections  II,  III,  IV, 
and  V  of  this  report. 
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a. 


0  <  a<  — 


The  energy  balance  for  T2  is 
^cond  j^2  +  ^cond  3^2  1  ^stored 


1  —  2 


3—2 


or 


K-«)  (T'”'Ti)  ,V^  (Ti'T;) 

(t  2  -  T2) 


P..  c.  T 
a  a  a  A  t 


A2 


where 


A  4> 
1-2 


A  cj> 
3-2 


A2  =  4> 

Letting 


R  +VV2  -  - 

i  4J 

•E->t 


Ta  •  a\  ^  (Ta  -  a) 

+  ‘ 


12 


v 

12 


6. 


k  At 

d 


^a  ca  'a 


T  2  »  ~  ^  +  ^  T2  t 


and  solving  for  T2 


t  Vl.-V 

I  e 

(/  vr-'V‘ 


b.  — -  <  a  <  xa 


M  -M 


f*-:) 


The  energy  balance  for  T2  is 
owl  j‘^2  ^cond  3l\  ^stored 


(M2) 


(343) 


(344) 

(343) 

(346) 


(347) 


1  —  2 


3—2 


(348) 


10-1 


Ra 

fT  n.;\ 

(Ta  '  a) 

l  1  m  x  2J 

(  wa  \  (r;  -  tJ 

Pa  cal  —  -  a  I  At 


r>4<)) 


whore 


.  ,  /„  Ta  “  a  \  (Ta  ;t) 


A  4> 

3-  2 


'(« £ 


Ta  \2  V 


72  +  T  +  TT 


a2  +  [(r +]>}*- 

Lotting 

ka  At 

a  ‘  Pa  ca  Ta2 
1 

and  solving  for  T2 


Ta  -  M  (3Ta  -  2a) 
4  /  4  48 


Z  =  R+£t2  , 


(5  50) 

(351) 

(352) 


A, •  »\  (?  -t)  -tt _ 

\  ’4  /  /  liiil)1  <tT-»  - 

_ y  <  /  * _ <#  . 


(~r  •  *)f*  •  ■*) 


T  ri  <r''r)^  v.  a.-a  T&-*. 

1  m  .  >  i  *  Tj  (lj  [  ■  ■  ,  *  T|  ■ 

L  '« •  -■■)  t»’a  •  \  ’»  /  /  T«  •  \  0’..  • 

Lr  ~  4  ’  18  _ tv  I  /  *  48 _ _ 

I,  I".-*-!1  ’•  /  i  ■  --  i 

V  I  /  48  V  4  /  ‘  48 


(353) 


Equations  (136)  and  (347)  are  similar  and  Equations  (139)  and  (353) 
are  similar.  If  a  multiplication  factor  is  substituted  into  Equation  (136), 
one  obtains  Equation  (  347).  The  same  is  true  for  Equation  (139)  to 
obtain  Equation  (353).  The  resulting  substitutions  and  equations  are: 


z  ■-  R  -  (XQ)^t2  . 


.  ‘  a 

c.  0  £  a  <  — 


(354) 


/  f4  *  4  t  I'  ,  •  il 

V  ■  <*«>  —7—1  •  •  ~ 

,  V 

'  '  i; 

7  *  l' 

A  V  4 


(355) 
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d. 


<  a  <  t 


-  *  a 


*/  ,  i 

/  y  1 

i 

f'r  lyfii  **  | 

Ti 

Tn,  iS, 

\  :  1  *  n 

V  ,x0)  77  * 

VT  '  */ 

f, .  rnt  v: u) 

( 

y 

l  'J  -1' 

_ 

X  ‘  i  1  4s 

V  •  1X01  4 

/  •  i» 

j 

f 

—  « 

/"* 

/  |l'*  ■ 

[?.  ♦  (XOl  -  1  . 

♦ 

(7 

**  -  i*Y 

\~ 

) 

V  1  /  4s 

y 

(XQ)  ^  J  • 

48 

(356) 


Factor  XQ  takes  on  the  value  of  +1  for  centripetal  flow  and  -1  for 
centrifugal  flow.  For  centripetal  flow  the  radius  to  the  original  outer 
surface  is  used,  while  for  centrifugal  flow  the  radius  to  the  original 
inner  surface  is  used.  Only  one  set  of  equations  is  compared  here; 
however,  all  of  the  equations  for  the  spherical  flow  can  be  arranged  in 
this  form. 


Due  to  the  nondimensional  terms  used,  it  is  also  possible  to  use 
a  -R  value  for  R  and  achieve  the  same  results.  Thus,  the  centripetal 
equations  presented  in  the  main  body  of  this  report  for  spheres  can  be 
used  in  their  present  form  to  calculate  centrifugal  heat  flow  effects  by 
inserting  the  radius  to  the  initial  inner  surface  as  a  negative  value. 


/ 
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Appendix  B 


COMPARISON  OF  THEORETICAL  TEMPERATURE  RESULTS  USING 
FINITE-DIFFERENCE  TECHNIQUES  FOR  FLAT  PLATES  DURING  ABLATION 


The  designer  or  analyst  always  encounters  a  basic  question  when 
using  numerical  techniques;  "How  well  do  the  results  from  these 
techniques  compare  with  data  obtained  from  exact  solutions?"  There 
is  only  one  exact  solution  that  may  be  readily  used  to  obtain  data  for 
comparison  with  data  generated  from  the  forward-backward  finite- 
difference  equations  described  in  this  report.  This  exact  solution  is 
for  a  semi- infinite  solid,  ablating  at  a  constant  rate  and  surface  tem¬ 
perature,  with  the  ablated  material  removed  from  the  surface  and  swept 
downstream.  This  solution  also  assumes  constant  thermal  properties. 


3 

The  exact  solution  is 

rp  rp  3.  X 

1  x  ~  1  CO  -  


where 


(357) 


Tx  =  Temperature  at  distance  x  in  from  the  exposed  surface. 
Tm  =  Melting,  ablating,  or  subliming  temperature. 

Too  =  Temperature  at  x  =  ®  from  the  exposed  surface, 
a  =  Ablation  rate. 

x  =  Distance  in  from  the  exposed  surface, 
a  =  Thermal  diffusivity  of  the  material. 


The  General  Thick  equations  were  used  to  obtain  temperature  data 
for  comparison  with  results  from  exact  solutions.  The  input  values 
for  the  exact  and  finite-difference  methods  were  the  following: 

ablation  rates,  a,  =0.1,  0.25,  0.4,  and  0.  5  mm/sec. 

Tm  -  2000°K . 

Tw  =  300°  K. 

Ta  =  0.001  m  =  1.0  mm. 

At  =1.0  sec. 


3H.  S.  Carslaw  and  J.  C.  Jaeger,  CONDUCTION  OF  HEAT  IN 
SOLIDS,  Second  Edition,  New  York,  New  York,  Oxford  University 
Press,  1959. 
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c  _ 

cl 


Pa 


a 


0.  4  kcal/Kg  -  ‘  K 
0.  56  kcal/m-hr-  K 


0.  25 


1000  Kg/m3 


9  x  10* 1  m2/hr 


The  steady- state  temperature  gradient  for  the  finite-difference 
method  was  obtained  by  raising  the  surface  temperature  of  a  semi¬ 
infinite  slab  to  the  ablating  temperature  at  time  zero  and  calculating 
(ho  surface  recession  and  temperature  distributions  for  a  sufficient 
tilin'  to  obtain  a  steady,  nonchanging  temperature  profile  (with  respect 
to  distance  from  the  receding  surface)  in  the  slab.  Comparisons  of 
these  steady- state  temperature  profiles  with  those  obtained  from  exact 
solutions  are  presented  in  Figures  52  through  55  for  four  different 
recession  rates.  The  temperature  data  comparisons  presented  in  the\se 
figures  show  that  the  accuracy  of  the  finite-difference  method  is  best  \ 
wln'ii  the  amount  of  material  removed  during  a  calculation  time  incre-  \ 
nient  is  small  in  relation  co  the  selected  incremental  distance  between  \ 
temperature  nodes.  For  exanip'e  with  the  |5  and  ablation  rate  held 
constant  at  0.25  and  0.  5  mm/sec  respectively,  the  accuracy  of  the 
finite-difference  method  is  improved  by  decreasing  the  At  and  t  as 
shown  in  Figure  55.  This  new  selection  of  parameters  reduces  the 
r  o  of  a  At/T  thereby  improving  the  accuracy  of  the  numerical 
approximations. 

Rased  on  the  temperature  gradient  comparisons  in  Figures  52 
through  55,  a  parameter  may  be  established  as  a  guide  in  selecting 
proper  inputs  which  will  result  i  acceptable  finite-difference  accuracies. 
If  the  finite-difference  temperatu  deviations  shown  in  Figures  52  and 
5  5  are  acceptable  and  those  sh  iwn  in  Figure  54  and  solution  No.  1  of 
Figure  55  „re  not  acceptable,  or  instance',  an  upper  limit  of  0.25  for 
a  At/T  is  established.  i\iathematically  this  criterion  is  stated  as 


a  At 


0. 


25 


(  558) 


r 


a  is  the  maximum  expec  ted  re  cession  rate-. 


Equation  (358)  can  be  written  in  another  form  by  considering  that 

PkAt  ,  k  .  . 

=  - =•  and  a  =  — .  This  equation  is 

pc T2  pc 

— P  <_  0.  25  (359) 

Equation  (359)  contains  the  term  aT/'a  which  is  equivalent  to  ax/a. 
in  Equation  (357).  If  this  term  is  too  large  the  relative  temperature 
difference  between  Tm  and  Tx  _  T  is  quite  large.  Based  on  the  con¬ 
ditions  considered  in  the  steady- state  temperature  comparisons  of 
Figures  32  through  35,  an  acceptable  upper  limit  for  the  finite- 
difference  criterion  aT/a  is  unity.  With  ax/a  1  the  temperature 
difference  between  Tm  and  T2  (located  one  T  from  Tm)  is  a  maximum 
of  approximately  63  percent  of  the  difference  between  Tm  and  the 
initial  equilibrium  temperature  of  the  slab. 

Satisfying  the  conditions  of  P  <_0.  5,  aAt/T  <_  0.  25,  and  aT/a  1 
before  performing  a  finite-difference  recession-condition  analysis 
aids  in  insuring  that  the  computed  heat  transfer  in  the  slab  will  be 
reasonably  accurate. 
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13.  Abstract  (Concluded) 

The  cylindrical  and  spherical  equations  are  presented  for  centripetal 
heat  flow  and  surface  recession.  Two  simple  methods  of  converting  the 
centripetal  equations  to  the  centrifugal  form  for  applications  to  structures 
such  as  blast  tubes,  rocket  motor  combustion  chambers,  and  nozzles  are 
discussed.  These  two  methods  involve  making  a  minor  number  of  sign 
changes  in  the  centripetal  heat-flow  equations. 

Attractive  features  of  the  ablation- conduction  method  described  in 
this  report  are  the  negligible  increase  in  required  computer  time  over  a 
nonreceding  case  when  all  other  parameters  are  identical.  Secondly, 
the  nonshifting  temperature  grid  prevents  confusion  in  interpreting 
computer  results  and  readily  lends  itself  to  automatic  plotting  techniques. 
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